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An important and well established area of quantum optics is the theory of Markovian stochas- 
tic Schrodinger equations (or by another name quantum trajectory theory). Recently stochastic 
Schrodinger equations have been developed for non-Markovian systems. In this paper we extend 
the current known stochastic Schrodinger equations for non-Markovian systems to include the po- 
sition unraveling. We also discuss and illustrate that this stochastic Schrodinger equation can have 
an interpretation under both the orthodox and the de Broglie-Bohm hidden variable interpretation 
of quantum mechanics. We conclude that only the de Broglie-Bohm hidden variable theory provides 
a continuous-in-time interpretation of the non-Markovian stochastic Schrodinger equation. 

PACS numbers: 03.65.Yz, 42.50.Lc, 03.65.Ta 



I. INTRODUCTION 

In nature we are more likely to find a system interact- 
ing with an environment than isolated. Hence a common 
problem in physics (especially in the area of quantum 
optics) is to model open quantum systems [1]. These 
systems consist of a small system (the system of inter- 
est) immersed in a lager system, which we refer to as the 
environment or bath. Due to the large Hilbert space of 
the bath it is convenient to describe the system by its 
reduced state. The reduced state is defined as 

Pred(t) =Trbath[|*W)(*WI], (1.1) 

where ^(t)) is the combined system-bath state, found 
from the Schrodinger equation. 

It has been shown [2, 3] by a projection-operator 
method that we can write a general (non-Markovian) 
master equation for the reduced state as 

i - /•* - 

dtpied(t) = --[H int (t),p Ied (t)}+ / fC(t,s)pred(s)ds, 
11 J to 

(1.2) 

where H- ln t{t) is a system operator in some interaction 
picture and Kit, s) is the 'memory time' superoperator. 
It operators on the reduced state p re d(t) and represents 
how the bath affects the system. The problem with this 
equation is that in general the effect of K{t, s) on p re d{t) 
can not be explicitly evaluated. 

Recently non-Markovian stochastic Schrodinger equa- 
tions (SSEs) [4, 5, 6, 7, 8, 9, 10] have been proposed which 
allow an alternative procedure for solving the reduced 
state. A non-Markovian SSEs is a stochastic equation 
for the system state \tfj z (t)), conditioned on some noise 
function z(t, s). (The double time argument in z will be 
explained in Sec. III.) The SSE has the property that 
when the outer product of \ifj z {t)} is averaged over all the 
possible z(t, s) one obtains p re d(t). That is, 

p ied (t) = E[\ip z (t))(ip z (t)\} > (1.3) 



where E[...] denotes an ensemble average over all possible 
z(t, s)'s. 

When using non-Markovian SSEs to solve the reduced 
state it turns out that in general we can not explicitly 
evaluate |f/> z (t))[7]. That is, we can only explicitly write a 
non-Markovian SSE for situations where a reduced state 
can be found exactly be other means. However as shown 
in Ref. [11] one can perform a post-Markov perturba- 
tion to the non-Markovian SSE, which allows approxi- 
mate results for p rc d(t) when the system is close to being 
Markovian. In Ref. [7] we illustrated a different per- 
turbation method that allows a perturbative result even 
when the system is strongly non-Markovian. However it 
is only valid for certain environment correlation functions 
(memory functions). 

Here we are not interested in how to solve current 
non-Markovian SSEs. The aim of this paper is to firstly 
extend the known non-Markovian unravclings (different 
functional forms of z{t, s)) to include what we label 
the 'position' unraveling, and then secondly to outline 
a non-local hidden variable interpretation of this non- 
Markovian SSE. This hidden variable interpretation is 
similar to the de Broglie-Bohm interpretation of quan- 
tum mechanics [13, 14, 15]. In other work we will supply 
a general hidden variable interpretation to include the 
coherent and quadrature unravelings [16]. These unrav- 
clings arc important to the quantum optics community 
as they have as their Markovian limit heterodyne and 
homodyne detection. The position unraveling which is 
presented here does not have a well defined Markovian 
limit. 



II. GENERAL DYNAMICS FOR 
NON-MARKOVIAN SSES 

The aim of this section is to outline the model used to 
develop non-Markovian SSEs. The results of this paper 
are only applicable to situations when the dynamics of 
the open quantum system can be described by the total 
Hamiltonian 
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H sys {t)®i + l®H hath + V, 



(2.1) 
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where H sys (t) = Hq + H(t) is the system Hamiltonian, 
and the bath is modelled by a collection of harmonic 
oscillators (for optical open quantum systems this corre- 
sponds to the electromagnetic field). In terms of dimen- 
sionless position (X k ) and momentum (Y k ) operators, the 
Hamiltonian for the bath is 



bath 



(2.2) 



where K is the total number of modes in the bath. These 
dimcnsionlcss operators have the commutator [Xj, Y k ] = 
iSj,k- 

The interaction Hamiltonian, V we assume is linear. 
By this we mean it has the form 



v = in 



EN 



(X k - iY k ) 
V2 



9k- 



(X k + i%) 



V2 



(2.3) 



where g k is the coupling strength of the k th mode to the 
system. 

For calculation purposes we define an interaction frame 
such that the fast dynamics placed on the state by the 
Hamiltonians Hq and -ffbath is moved to the operators. 
The unitary operator for this transformation is 

U (t,to) = e ^ i ^ u ^ + ^®^ bath ^ t ~ t "^ h . (2.4) 
This allows us to write the Schrodingcr equation as 



d t m)) 



h 



[^intW+Ant(t)]|*W), 



(2.5) 



where H m t{t) is H(t) in the interaction frame and 

K 

V int (t) = [Lg* k e iQ * 



i k (t-t ) ( X k-iYk) 



k 

-L^g k e 



V2 

-m k (t-t a ) + iYk) 



V2 



(2.6) 



with Q k = cu k — CI. Here we have finally restricted the 
form of Hq to be such that L in the interaction frame 
simply rotates in the complex plane at frequency f2. That 
is L int (t) = L e -«(*-*o), 



A. Quantum Measurement Theory 

In open quantum systems a measurement is usually 
performed on the bath rather than directly on the system. 
Due to the entanglement between the bath and the sys- 
tem the measurement on the bath results in an indirect 
measurement of the system. For the position unravel- 
ing the measurement performed (on the bath) is to mea- 
sure the set of dimensionless position operators {X k }. 
To mathematically describe this measurement process we 
consider the if-observables 

{X k } = {({x k },it {xk} = \{x k })({x k }\ ® l sys )}. (3.1) 

Here the eigenvalues {x k } represent the actual results of 
the measurement and the eigenstate is the multi- 

mode state the bath is projected into upon measurement. 
Because this theory is done in the interaction frame the 
measurement is not exactly a position measurement. It is 
position defined in the interaction (rotating) frame which 
in the Schrodingcr picture (stationary frame) represents 
a measurement which cycles between position and mo- 
mentum. 

To find the state of the bath and the system after the 
measurement we decompose the projector into measure- 
ment operators, M{ Xk y, by 7r {xk} = M f {xk} M {xk} . With 
this measurement operator the combined state after a 
measurement at time t, which yielded results {x k } is 



(3.2) 



where P({x k },t) is the probability density and is defined 
by 

P({x k },t) = (*m {xk} ® i sys |*(i)>. (3.3) 

For the rank-one projector, |{xfc})({a;fc}|, the measure- 
ment operator which satisfies the above decomposition is 
Mi Xk \ = \{n k })({xk}\, where is arbitrary and is 

the state the bath is left in after the measurement. For 
most optical situations this will be a vacuum state |{0fe}) 
(as in optical measurement the detector usually absorbs 
the photon). This results in the combined state becoming 
\^{x k }(t)) = |{0fc})|Vw(*)>, wher e |^W(*)> is denne d 



III. DERIVING THE POSITION 
NON-MARKOVIAN SSE 

Under the orthodox view [12] of quantum mechan- 
ics the quantum state upon measurement undergoes a 
change which is consistent with the measurement re- 
sults. This change in state has been termed a 'collapse'. 
Whether this collapse is a real physical event or repre- 
sents an update in the observer's knowledge will not be 
consider here. The central point is that under this view 
our observation causes the wavefunction for the state to 
collapse. The standard way to model this collapse is to 
use quantum measurement theory (QMT) [17]. 



(MM*)) 
VPWkUY 



(3.4) 



which we label a conditioned system state. This condi- 
tioned system state has the property that 



Pred(*) = J P({x k },t)\ij {xk} (t))(lP {xk} (t)\d{x k }. (3.5) 

If we consider {x k } to be a set of random variables 
{x k (t)} chosen from the distribution P({x k },t) then we 
can rewrite the reduced state as 

Prcd (t) = £[|Vw*)}(*)W{**(*)}(*)l]> ( 3 - 6 ) 
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where E denotes an average over P({xk},t). This is the 
requirement outlined in Sec. I for a solution to a SSE. 
This suggests that the time derivative of \tp{x k (t)}(t)) will 
be a non-Markovian SSE. To calculate this we need to be 
able to generate a self consistent differential equation for 
I V'fxfcft) }(£))• This is complicated and to simplify this pro- 
cedure we introduce linear quantum measurement theory 
(LQMT). 



B. Linear Quantum Measurement Theory 

LQMT uses the same principles as QMT except we 
use an ostensible distribution, A({xk}), in place of the 
actual distribution [18, 19]. Using this distribution the 
linear conditional state is 



as 



fe fe }(*)> 



(3.7) 



The bar above this linear state signifies that the state is 
unnormaliscd. As before the reduced state can be written 



PrcdW = / A({x k })\$ {Xk} (t))$ {xk} (t)\d{x k }. (3.8) 



Thus p Ied (t) = E[\ip {Xk{t)} {t))(ip {Xk{t)} {t)\} where E de- 
notes {x k (t)} is chosen from the ostensible distribution 
A({xfc}). Because A({xk}) is time independent {xk(t)} = 
{x k (t )} for all t and d t \4'{ Xk }(t)) = d t \i>{ Xk }(t)). 

Using Eq. (3.7) the time derivative of |"0{ a;) .}(i)) be- 
comes 



d t |Vw(i)) 



({x k }\dt\*(t)) 

VMM) 



(3.9) 



Substituting into this Eq. (2.5) and the Hamiltonian de- 
fined in Eq. (2.6) the differential equation for the linear 
system state becomes 



J 



dt\i>{ Xk} (t)) = [-j.H int (t)+ yk 



m k (t-t ) 



V2 

^2 k (Lg* k e ia ^-^ +frg k e- in *V-^) d Xk {{x k }\*(t)) 



IVW(*)) 



V2 



Choosing the ostensible distribution to be 



A(W) = |({.T fc }|{0 fe })| 2 = W'-^jf 1 : 



(3.10) 



(3.11) 



results in 

dM {Xk} {t)) 



iQ k (t-t ) 



9x k 



lVwW>- 



(3.12) 



r 



Note this ostensible distribution was chosen for simplicity 
(it is equal to the real distribution at t = to). 

To define a SSE from Eq. (3.12) we make the random 
variables substitution x k — > x k (t). This transforms the 
linear state as |^{ a:fc }(i)) — * |^{x fe (t)}(A))- F° r the linear 
SSE the random variable x k (t) must satisfy E[x k (t)\ = 
and E[x k (t)xi(t)] = 6 k ,i/2 for all t (due to Eq. (3.11)). It 
can be still argued that, even with this substitution, this 
is not strictly a SSE as to solve this we must simulta- 
neously solve for all possible {x k (t)} (due to the partial 



derivative in Eq. (3.12)). To make it a SSE we need to re- 
place the derivatives by operators, but before we do this, 
to keep with past literature, it is convenient to define a 
noise function z{t, s) as 



z(t,s) = y^^g k V2x k (t)e~ 



iQk(s-to) 



(3.13) 



With this definition we can define the following two bath- 
correlation function 



E[z(t,s)z*(t,s')} = ^g h gt2e- ia ^'-^ +ia ^-^E[x k (t)xi(t)] 
k,l 



(3.14) 



4 



J2\9k\ 2 e- m ^-'^=a(s-s'), 



E[z(t,s)z(t,s')\ = J2gle- m ^ +s '- 2 ^= 7 (s + s') 
k 



and the K partial derivatives become 



Sz(t,s) 



5z*(t,s) 



where 5 z »(t jS ) and 5 z *u tS \ are functional derivatives. In terms of the noise function, Eq. (3.12) becomes 



d t \i> z {t)) = {-^H int (t)+Lz*(t,t)-L £ [a*(t-s)5 z(t , s) + 7* (* + 
+7(t + s)<J,( tl »)]<k}|&(*)>- 



1 r 



(3.15) 
(3.16) 



(3.17) 



(3.18) 



Assuming the initial combined state to be |\I>(io)) = 
|{0}fe)|V'(^o)) it is easily shown that \tp z {t )) = |V>(*o)> 
thus the functional derivative with respect to z(t, s) in 
the above equation will always have zero contribution. 
That is, we can rewrite Eq. (3.18) as 



d t \Mt)) = {-iH int ( t )+Lz*(t,t)- [i 7 *(t + s ) 
+Pa{t- s)]5 z * (ttS) ds}\i> z (t)). (3.19) 



This equation appears nicer than Eq. (3.12), but it is es- 
sentially the same equation and we still have the problem 
of representing the functional derivative by operators. To 
do this we make the two ansatzen, 



A z (t)$ z (t)) = / dsa(t-s)5 z * (t!s) \ip z (t)), (3.20) 



B z {t)Hz{t)) 



dsj*(t + s)6 z * (t)S) \j> z (t)). (3.21) 



C. The Non-Markovian SSE 

In the above section we have presented the linear SSE 
for the position unraveling. To generate the actual non- 
Markovian SSE from this linear equation we have to first 
of all generate the set of random variable {x k (i)}, which 
obey the real distribution P({x k },t). These in general 
will be time dependent. To do this we use a Girsanov 
transformation [4, 20] to relate the real probability dis- 
tribution with the ostensible one, that is 

P({x k },t) = A({x k })(iJ {xh} (t)\$ {xk} (t)). (3.23) 

Taking the time derivative of this and using Eq. (3.12) 
gives (after some manipulating) 



d t p{{ Xk \,t) = -$>»>{^ {Xfc} (t)|[£s 



* ififc(i-io) 



fee 



-ifi fc (i-i ) 



]|VwW> 



'P({xk},t)} 



(3.24) 



With these ansatzen the linear non-Markovian SSE for 
the position unraveling is 

d t \Mt)) = {- l -H int {t)+Lz*{t,t)-LB z {t) 

-DA z {t)]\i) z (t)). (3.22) 

In general these two ansatzen can not be solved, but as in 
coherent and quadrature unravelings [7, 11] there are per- 
turbation techniques which could be extended to include distribution (t )\Tt Xk \^ (t )) , which for an initial com- 



This is effectively a drift equation for the probability den- 
sity. It has associated with it the following set of differ- 
ential equations 

dkx k {t) = [(L} t g* k e m ^-^ + (L^ t g k e- in ^-^}/V2 

(3.25) 

where (L) t = (■ip {xk(t)} (t)\L\if; {xk{t)} (t)) . Integrating this 
gives 

x k (t) = x k (t Q ) + f dt'[(L) t ,gie in ^'-^ 

(3.26) 

where x k (to) is the random variable associated with the 



+ (L^g k e~ m ^'~^]/V2, 



this unraveling and would allow a pcrturbative solution. 
This we leave for later work. 



bined state of the form \^(t )) = \{O k })\ip(to)) is equiv- 
alent to the above ostensible distribution, Eq. (3.11). 
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The noise function for the real distribution, defined in 
Eq. (3.13), becomes 

z(t,s)=z(t ,s) + f dt'{(L} t ,a(s-t') + (fr) tn (s + t')] 

(3.27) 

where z(to,s) is equivalent to the noise function used in 
the linear non-Markovian SSE. 

To generate the non-Markovian SSE we use Eqs. (3.4) 
and (3.7) to rewrite the conditioned state \ip{x k (t)}(t)) = 

(*)) !{**=**(*)> as 

J {^{x k (t)}{t)\^{x h (t)}{t)) 



As long as {xk{t)} in \^{x k (t)}{t)} are the random vari- 
ables associated with the real distribution, Eq. (3.26), 
this will obey the correct statistics defined in Eq. (3.6). 
Using Eq. (3.22) and the standard definition for a to- 
tal time derivative gives the following SSE (after some 
manipulation) 



dt\M*)) = {-^H int (t) + {L-(L) t )z*(t)-(L- 
- (P) t )A,(t))t}\Mt)h 

where \ip z (t)) = |-0{x fe (t)}(O)- Here we have again used 
the above ansatzen to remove the partial derivatives with 
respect to {xk}- Under the orthodox interpretation, the 
solution of this equation at time t is the state of the 
system given that a measurement has been performed on 
the bath at that time and yielded results {xk(t)}. Since 
this measurement would change the bath state, the future 
evolution of the system would not be the same as if the 
measurement had not been performed. In other words, 
the solutions at different times correspond to different 
physical situation. Thus the linking of these solutions to 
make a trajectory for the system state is a fiction. 

IV. SIMPLE APPLICATION 

To illustrate that this is a correct unraveling for a non- 
Markovian open quantum system, in this section we ap- 
ply this theory to a simple model. The model is a two 
level atom (TLA) coupled linearly to a single mode bath 
(K = 1) with no detuning. To get an exact solution 
for comparison we first of all solve the total Schrodingcr 
equation for the combined state. For this system the 



(L) t )B z {t) + ((L - (L)t)B z (t)) t - {p - P) t )A z {t) 

(3.29) 

I 

Schrodingcr equation in the interaction frame is 

\^t)) = {g*aa) -g^a)m)) (4.1) 

where L = a = \b)(e\ is the lowering operator for the 
TLA and |e) and \b) are the excited and ground states of 
the TLA. The combined state can be written in terms of 
photon number states and atomic states as 

!*(*)>= E E c *»n- ( 4 - 2 ) 

s=e,b n 

Substituting this into Eq. (4.1) and using the fact that 
since the bath state is initially in a vacuum state the 
only nonzero amplitudes are c e ,0) Q>,o and cj,,i, we get 
the following solutions (with to = 0) 

c e ,o(i) = ce,o(0)cos[|ff|(t-td)], (4.3) 
c 6> o(i) = c M (0), (4.4) 
c b ,i(t) = c e ,o(0)sin[| 5 |(i-i )]e- ie , (4.5) 

where 8 is the argument of the complex coupling constant 
g. Thus the reduced state is simply 



Pred(t) = < o (i ) cos 2 [| ff |(t - *o)]|e><e| + {cg i0 (t ) + < (*o) sin 2 [j 5 |(t - to)]}\b)(b\ 

+ [c e , (£o)< (t )|e}<&| + c* O (to)ct,o(*o)|fc)<e|] cos[\g\(t - to)]. (4.6) 

The general non-Markovian SSE defined in Eq. (3.29) for this simple system becomes 

dt\1>,(t)) = {(* - (&)t>*® ~ (* - (*)t)B z (t) + ((& - (°)t)B z {t))t - (^ - (^)t)A z (t) 

+ ((v f -(^)t)A z (t)) t }\Mt)), (4.7) 
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where 



z(t) = gV2xi(t) 



(4.8) 



and 



xx (t) = Xl (t ) + / dt'(g*(a) t , + g(^) t ,)/V2. (4.9) 



The value of the random variable xi(to) is determined by 
the initial distribution 

P(x 1 ,t )=A(x 1 ) = eXp{ - x2l) . (4.10) 



To find functional form of the two operators A z (t) and 
B z {t) we use Eqs. (3.20) and (3.21) and assume 

*«•(.)!&(')> = f(t,s)a\Mt)), (4.11) 

This results in A z (t) = A z (t)a, and B z (t) = B z (t)a, 
where 



A z {t) = / ds\g\ 2 f(t,s), (4.12) 

-'to 
ft 



dsg* f(t,s). 



(4.13) 



That is A z (t) and B z (t) are independent of the noise 
function for this example. Taking the time derivative of 
these equations gives 



AM = \g\ 2 + / \g\ 2 f(t,s)d 



(4.14) 



B z (t) = g* 2 + / g* 2 f(t, S )ds. (4.15) 

Jto 

Here we have used the fact that f(t,t) = 1 [6]. To find 
the form of f(t,s) we firstly need to define the linear 
non-Markovian SSE for this system. It is 

d t \Mt)) = [&z*{t) - A z (t)^a}\Mt)), (4.16) 

where z(t) = g^/2x\(t) with x\(t) in the linear case being 
defined to have a value given by A(xi) [see Eq. (4.10)]. 
Secondly we use the following consistency conditions 

S zH s)d t \Mt)) = (4.17) 

Substituting Eq. (4.16) and Eq. (4.11) into this equation 
gives the following differential equation 

f(t,s) = f(t,s)A z (t), (4.18) 

which when substituted into Eqs. (4.14) and (4.15) gives 

A z (t) = \g\ 2 +Al(t), (4.19) 




B z {t) = g* 2 + B z {t)A z {t). 



(4.20) 



FIG. 1: An example trajectory for the position non- 
Markovian SSE for a TLA in a single mode bath. Also shown 
is the real (solid) and imaginary (dotted) part of the noise 
function. Note the imaginary part for this system is equal 
to zero for all time. All calculations where done numerically 
with g = 1, a time step size of 0.0001 and an excited state 
initial condition. 



Solving this set of coupled differential equation with the 
initial conditions A z (to) = B z (to) = gives A z (t) = 
|ff|tan[|ff|(t-to)], and B z (t) = e - ae \g\t&ii[\g\(t - t )}. 
Thus Eq. (4.7) is now numerically solvable. An exam- 
ple solution for \ip z (t)) is shown in Fig 1. Here it is 
observed that unlike a Markov trajectory the evolution 
is smooth. This is because with only one bath mode 
the bath correlation time is non-zero. In this figure the 
quantum state is displayed in the Block representation, 
that is x(t) = (a)t + (cf^)t, y(t) = — i(&)t + i{^)t and 
z{t) = ((T^<j)t — (cfa^)f To show that this equation does 
reproduce the reduced state, the difference between the 
ensemble average of 1000 trajectories and the reduced 
state, Eq. (4.6), is shown in Fig. 2. It is observed that 
the ensemble averages does reproduce the reduced state 
within statistical error. 



V. A HIDDEN VARIABLE INTERPRETATION 
OF THE POSITION NON-MARKOVIAN SSE 

In this section we present an alternative interpretation 
of the position non-Markovian SSE. The hidden variable 
interpretation we use is similar to the de Broglie-Bohm 
interpretation. The only difference is our position vari- 
ables are 'positions' in a reference frame which rotates in 
phase space (the interaction frame). This results in our 
trajectories for the positions to be different from the stan- 
dard de Broglie-Bohm trajectories [13, 14, 15]. Moving 
to this frame simply means that the position states as we 
have defined them actual correspond to, in a stationary 
frame, a state which rotates between position and mo- 
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FIG. 2: The difference between the ensemble average of 1000 
trajectories for the position unraveling and the master equa- 
tion for a TLA in a single mode bath. Other details are as in 
figure 1. 



mentum. This is not a problem provided the Schrodingcr 
equation used to develop the trajectories is in the inter- 
action (rotating) frame. 

In the standard dc Broglic-Bohm theory, position has 
a reality prior to the measurement and its trajectory is 
deterministic assuming its initial position is known. To 
account for the weirdness of quantum mechanics its tra- 
jectory can be very non-classical. Specifically, the tra- 
jectory is influence by a extra potential, the quantum 
potential [14, 15]. This quantum potential depends on 
the solution of the Schrodingcr equation, and is why the 
wavefunction is called the guiding wave [13, 14, 15]. Thus 
under this interpretation the wavefunction is a real entity. 

In this paper we will not be introducing the quantum 
potential as we can describe the trajectories with refer- 
ence to only the wavefunction as a guiding wave. To 
mathematically describe this interpretation we use the 
fact that since the probability density is continuous in 
time and configuration space {x k e $l K } and is a con- 
served quantity, it must obey a continuity equation, 



d t P({x j },t) = -^2d x j k ({ Xj },t), 



(5.1) 



where jk({xj},t) is the current density and is related to 
the velocity field by jk({xj} 7 t) = P({xj}, t)v k ({xj} 7 1). 
Each individual trajectory is then found by 



d t x k (t) = V k ({xj},t)\{ Xj=XjW y. 



(5.2) 



Thus to find this trajectory we need to know either 
Vk({xj},t) or P({xj},t) and jk({xj},t). In previous lit- 
erature the method is normally to find P({xj},t) and 
jk({xj}, t). Here, however, we propose that v k ({xj}, t) is 



given by 



Vfc({x 3 -},t) 



Ite[(*(*)|{^})({^}|v fc (t)|*W>] 



(5.3) 



(*(*)lte}><{*i}|*(*)> 

where v k (t) is called the velocity operator. It is given by 



v fc (t) = -~[X fe ,Vi nt (t)]. 



(5.4) 



Note this is only valid for Hamiltonians with terms which 
are at most quadratic in Y (this is proven in Ref. [16]). 
This is not a problem as in nature all fundamental Hamil- 
tonians are of this form. 

Up until now we have not said anything about com- 
bined systems. If we now include an extra system, the 
system of interest, but only calculate trajectories for the 
bath positions, then the following differs slightly from 
the conventional de Broglie-Bohm theory Defining the 
system state as 



<{*,}!*(*)} 



the velocity field becomes 



(5.5) 



v k ({ Xj },t) =Re{(^ {Xj} (t)\v k ({x J },t)\4> {xj} (t))], (5.6) 
where 



v fc({Zj}, Wfa^IW) = ;)|T|/J . |f ffff ff= 



V(*(t)|{»i})({x 3 -}|*(t)) 



(5.7) 



Note the arrow defines the direction of operation 
[vfcdccj}, t) will contain partial derivatives the act in the 
direction of the arrow]. Substituting the velocity field 
[Eq. (5.6)] into Eq. (5.2) results in following differential 
equation for the bath positions 



d t x k (t) = Re[(tp {x (t)}{t)\v k {{xj (*)}, t)\ip {x ( t )} (*))], 

(5.8) 

where \4>{ Xj (t)} (*)} = IV^} (t))\ {Xj=Xj ( t)} . 

We now make the link that Eq. (5.5) is mathematically 
equivalent to Eq. (3.4) even though they have different 
physical interpretations. Since in Sec. Ill we showed 
that the position non-Markovian SSE is derived from this 
equation, then the time derivative of Eq. (5.5) must also 
give the correct non-Markovian SSE. Thus under this hid- 
den variable interpretation the solution of the position 
non-Markovian SSE is a real entity which exists for all 
time and guides the trajectories for the bath positions 
{xk(t)}. It gives the 'conditioned' system state at all 
times, as under the de Broglic-Bohm interpretation the 
bath positions exist even when we do not measure the 
bath. 

To show that Eq. (5.8) does give the same trajectories 
as used in the non-Markovian SSE derived from QMT, we 
used the above definition for Vfc(t) and apply the Hamil- 
tonians defined in Sec. II to them. This gives 



v k (t) = \g* k e in > t L + 



9ke 



(5.9) 
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which results in a velocity field of the form 

vfc(fo},i) = [^e^ 4 <V> { ^ } (i)|L|Vw(*)> 

+.g fce -^ t (^ j} W| J L t |Vfe}(t))]/V2- 

(5.10) 

Thus the actual trajectories are found from 

dtx k (t) = [{fytgie* 1 ** + (Lt) t5fce -^ 4 ]/V2 (5.11) 

where (L) t = (ifj {x . {t)} (t)\L\ip {x . {t)} (t)), which is the 
same as Eq. (3.25). 

VI. DISCUSSION AND CONCLUSION 

In this paper we have firstly extended the theory of 
non-Markovian SSE to include the position unraveling. 
In Sec. IV we applied this theory to a simple system, a 
TLA coupled linearly to a single mode, to show that the 
position non-Markovian SSE does correctly average to 
give the reduced state. Although this unraveling does not 
have a well defined Markovian limit, it is still interesting 
to look at as it allows us to investigate intcrpretational 
questions in quantum mechanics. This is because this 



unraveling can be interpreted under cither the orthodox 
interpretation or by a non-local hidden variable theory 
which is similar to the de Broglie-Bohm theory. 

Under the orthodox theory the solution of the non- 
Markovian SSE at time t is the state the system will be 
in if a measurement was performed on the bath at time 
t and yielded results {xk(t)}. Thus the non-Markovian 
SSE is simply a numerical tool for calculating the correct 
conditioned system state. In the dc Broglie-Bohm hidden 
variable interpretation the non-Markovian SSE is an evo- 
lution equation for the system state (a real state) which 
guides the trajectories for the baths individual positions 
{xk(t)}. Unlike the orthodox theory, these variables exist 
even when the bath is not measured. We conclude that 
if a non-trivial (non-numerical tool) interpretation is to 
be given to this non-Markovian SSE then we must con- 
sider the dc Broglie-Bohm hidden variable interpretation 
of quantum mechanics. 
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